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I.  ABSTRACT  (Maximum  200  words) 

Groundwater  contamination  by  substances  that  are  immiscible  with  water  (e.g.,  gasoline,  jet  fuels, 
chlorinated  solvents)  is  well-documented  as  a  commonpleu:e  occurrence.  The  release  of  such  nonaqueous 
phase  liquids  (NAPLs)  leads  to  a  multiphase  system  in  which  contaminants  can  exist  in,  move  within,  and 
be  transport  among  a  NAPL  phase,  an  aqueous  phase,  a  gas  phase,  and  solid  phase.  A  conunon  assumption 
until  now  is  that  each  chemical  species  is  in  local  equilibrium  among  all  the  phases  present  in  a  subsurface 
environment. 

This  research  project  investigated  the  propriety  of  the  local  equilibrium  assumption  through  conduct 
of  fundamental  laboratory  studies  and  mathematical  modeling  of  interphase  mass  transfer  in  multiphase 
porous  media  systems.  The  experimental  results  of  this  work  were  used  to  develop  interphase  mass  transfer 
rate  coefficient  correlations  for  NAPL-aqueous  and  aqueous- vapor  systems,  and  to  determine  correlations  of 
NAPL  morphology  at  residual  saturation.  The  modeling  portion  of  this  study  incorporated  the  experimental 
results  into  one-  and  two-dimensional  multiphase  flow  and  transport  models — in  order  to  investigate  expected 
deviations  from  equilibrium  conditions  at  the  field  scale.  Experimental  and  modeling  results  showed  that 
deviations  from  equilibrium  are  common  for  vapor-aqueous  phase  mass  transfer  and  may  be  important  for 
NAPL-aqueous  phase  mass  transfer,  especially  in  heterogeneous  systenvs. 
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Statement  of  Problem  Studied 


The  groundwater  resources  of  the  United  States  have  been  extensively  contaminated 
by  a  wide  variety  of  organic  pollutants.  This  widespread  contamination  has  led  to  a  need 
for  a  better  understanding  of  the  processes  that  impact  the  transport  and  eventual  fate  of 
such  contaminants.  A  better  understanding  of  these  fundamental  processes  is  necessary  to 
accomplish  more  accurate  assessment  of  health  risks,  and  more  economical  and  effective 
design  of  aquifer  renovation  measures. 

One  of  the  most  common,  and  most  challenging,  groundwater  contamination  problems 
results  from  the  spill  of  a  substance  that  is  immiscible  with  water.  Immiscible  contami¬ 
nants  (nonaqueous  phase  liquids  or  NAPLs)  are  frequently  introduced  into  ground vvaler 
systems  by  leaking  underground  storage  tanks,  the  associated  appurtenances,  or  other  di¬ 
rect  releases.  Once  released,  an  immiscible  solvent  begins  a  complex  path  of  transport  and 
transformation.  Accurate  description  of  the  resulting  multiphase  groiuidwater  system  is  a 
requisite  for  describing  the  time  dependent  mass,  and  distribution  of  contaminant  mass — 
an  important  objective.  To  this  end,  a  need  exists  to  better  imderstand  interphase  mass 
transfer  phenomena  in  such  systems. 

The  treatment  of  multiphase  systems  has  traditionally  been  within  the  realm  of  the 
oil  industry,  but  the  need  to  protect  groundwater  quality  has  changed  this  division.  The 
motivation  of  the  oil  industry  work  has  been  to  remove  as  much  of  the  immiscible  product 
from  the  subsurface  as  possible,  and  this  work  has  usually  neglected  the  interphase  mass 
transfer  of  components  of  the  immiscible  organic  phase  to  the  aqueous  phase.  The  study 
of  multiphase  groundwater  systems  is  oriented  differently — to  determine  the  impact  of  a 
NAPL  on  the  aqueous  phase. 

The  motivation  for  this  investigation  was  the  typical  treatment  of  solute  transport  in 
multiphase  systems  as  linear  equilibrium  processes.  Significant  evidence  existed  to  suggest 
that  the  assumption  of  a  rapidly-attained,  linear  equilibrium  between  phases  is  a  gross 
oversimplification  for  some  situations.  This  study  sought  to  quantitatively  elucidate  the 
effects  of  nonequilibrium  mass  transfer  conditions  in  subsurface  systems. 

The  goal  of  the  project  was  met  through  conduct  of  fundamental  laboratory  investi¬ 
gations  of  multipheise  mass  transfer  processes,  which  are  operative  in  subsurface  systems. 
This  investigation  focused  on  the  collection  and  analysis  of  a  detailed  set  of  quantitative 
data  on  the  final  phase  equilibrium  and  the  rate  of  approach  to  that  equilibrium  for  mul¬ 
tiphase  mass  transfer  reactions.  NAPL-aqueous,  NAPL-solid,  and  aqueous- vapor  phase 
processes  were  given  primary  consideration. 
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Summary  of  Most  Important  Results 


Overview 

As  noted  in  the  publication  section  of  this  report,  a  large  quantity  of  significant  and 
publishable  new  information  was  developed  during  this  research  project.  This  section  will 
briefly  summarize  the  most  important  major  findings  of  this  work.  Details  of  procedures 
used,  data  gathered,  and  full  discussions  of  data  analysis  and  interpretations  can  be  found 
by  reference  to  the  project  documents  produced  to  date  and  other  documents  currently 
in  preparation.  This  brief  summary  will  be  grouped  into:  NAPL-aqueous  phase  mass 
transfer,  residual  saturation  morphology,  and  aqueous-vapor  mass  transfer. 

NAPL-Aqueous  Mass  Transfer 

As  is  well  known  now,  once  an  immiscible  fluid  is  released  in  the  subsurface  it  tends 
to  move  vertically  tlirough  the  unsaturated  zone,  leaving  a  portion  of  the  pore  space 
filled  with  NAPL  that  is  held  inunobile  by  capillary  forces.  Thus  over  a  relatively  short 
time  period — compared  to  the  total  time  the  multiphase  fluid  persists  in  the  subsurface-  - 
the  multiphase  fluid  becomes  trapped  in  a  state  of  residual  saturation.  A  major  route 
of  subsurface  contamination  occurs  through  dissolution  of  the  trapped  NAPL  into  the 
aqueous  phase  that  flows  through  the  N  A  PL-contaminated  region. 

A  significant  portion  of  this  work  dealt  with  the  measurement  amd  modeling  of  the 
NAPL-aqueous  phase  dissolution  process.  This  is  significant  because  of  the  large  number 
of  sites  that  are  contaminated  with  fluids  that  are  immiscible  with  water.  A  paper  by 
Miller  et  al.  (1990)  details  the  experimental  procedures  used  to  measure  the  dissolution 
process,  the  development  of  dimensionless  variables  that  characterize  the  process,  and 
the  derivation  of  an  empirical  relationship  that  allows  for  the  estimation  of  mass  transfer 
characteristics  as  a  function  of  the  rele\’ant  physical  and  chemical  properties  of  a  system. 

Mayer  and  Miller  published  two  papers  on  the  development  and  application  of  a  gen¬ 
eral  multiphase  flow  and  transport  model.  The  NAPL-aqueous  phase  dissolution  process 
was  include  in  this  model  and  the  laboratory-derived  constitutive  relationship  was  used  to 
estimated  typical  mass  transfer  coefficients. 

The  results  of  both  the  experimental  and  modeling  work  show  that  NAPL-aqueous 
phase  mass  transfer  is  relatively  fast  compared  to  the  space  and  time  scales  of  concern  in 
typical  subsurface  systems.  Using  typical  parameters  for  groundwater  flow,  a  groundwater 
phase  moving  through  a  zone  of  uniformly  trapped  NAPL  would  be  expected  to  reach 
equilibrium  with  the  NAPL  phase  within  a  fraction  of  a  meter.  This  suggests  that  for  such 
conditions  equilibrium  between  the  water  and  the  NAPL  phase  is  a  reasonable  assumption. 

An  imanswered  and  potentially  very  important  question  coming  from  this  work  is  the 
effect  of  subsurface  heterogeneity  on  achievement  of  equilibrium.  That  is  heterogeneity 
of  both  the  NAPL  distribution  and  the  porous  media  system  would  be  expected  to  cause 
flow  by-passing  of  the  most  highly  contaminated  region.  This  could  lead  to  significant 
deviations  from  equilibrium  conditions.  Evaluation  of  this  point  will  require  additional 
future  research. 
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Residual  Saturation  Morphology 

As  was  noted  in  the  previous  section,  heterogeneity  is  expected  in  terms  of  the  distri¬ 
bution  of  trapped  NAPL  in  regions  at  residual  saturation.  This  natmal  variation  woiild 
be  amplified  by  variability  in  the  porous  media  characteristics.  The  morphology  of  the 
trapped  NAPL  (i.e.,  size  of  the  trapped  NAPL  regions,  volume  to  svurface  area  ratio,  and 
proximity  to  other  trapped  NAPL  regions)  will  bear  importantly  on  the  rate  at  which  the 
NAPL  will  be  dissolved  into  a  flowing  groundwater  phase. 

Experiments  were  performed  using  a  two-dimensional  micro-model  to  relate  physical 
and  chemical  factors  of  multiphase  systems  to  the  NAPL  distributions  noted  at  residual 
saturation.  These  studies  showed  that  the  viscosity  and  density  of  the  NAPL  were  very 
important  factors,  while  the  velocity  of  the  flowing  groundwater  had  a  minion  i’  effect 
of  the  NAPL  morphology — for  typical  ranges  of  groundwater  velocities.  This  suggests 
that  pump-and-treat  methods  will  not  be  effective  in  mobilizing  trapped  NAPLs  as  a  free 
phase.  These  results  also  suggest  that  as  the  density  of  a  NAPL  increases  the  surface 
area  to  volume  ratio  of  the  trapped  NAPL  will  increase,  thereby  favoring  local  equilibriiun 
conditions.  More  complete  experimental  and  interpretation  results  of  this  work  are  given 
in  a  manuscript  by  Mayer  and  Miller,  which  will  soon  appear  in  Transport  in  Porous 
Media. 

Important  questions  that  remain  to  be  answered  are:  the  effect  of  d'mensionality  of 
the  experimental  system  on  the  morphology  of  the  trapped  NAPL;  an  .  the  derivation 
of  the  characteristics  of  a  representative  elementary  volume,  which  is  required  to  perform 
meaningful  experimental  studies,  as  a  function  of  physical  and  chemical  system  parameters. 

Aqueous- Vapor  Mass  Transfer 

Mass  transfer  from  the  aqueous  phase  to  the  vapor  phase  is  important  not  just  in 
multipheise  systems  but  for  contaminant  fate  and  transport  for  any  volatile  organic  solute. 
This  issue  has  become  increasingly  important  in  recent  years  as  vacuum  extraction  has 
become  a  common  method  of  subsurface  restoration. 

A  series  of  experiments  were  performed  in  this  work  to  determine  the  rate  at  which 
mass  transfer  occurs  from  the  aqueous  phase  to  the  vapor  as  a  function  of  key  system 
variables.  The  results  of  this  work  were  interpreted  using  two  modeling  approaches:  a 
simple  analytical  model  that  ignored  longitudinal  dispersion,  and  a  two-dimensional  finite 
element  model  that  included  the  effect  of  longitudinal  dispersion. 

The  key  results  of  this  phase  of  the  work  were: 

1.  The  liquid-phase  mass  transfer  coefficient  of  toluene  in  water  (used  as  a  model  com- 
poimd)  increased  about  one  order  of  magnitude  for  every  two  order  of  magnitude 
increase  in  the  liquid-phase  velocity. 

2.  Liquid-phase  resistzmce  was  found  to  dominate  the  overall  mass  transfer  resistance  for 
toluene. 

3  Field-scale  simulations  showed  that  imder  ambient  conditions,  significant  portions  of  a 
subsxuface  system  would  not  be  described  accurately  using  an  equilibrium  assiunption. 

4.  Field-scale  simulations  of  a  vapor-extr2u;tion  scenario  showed  that  a  typical  aquifer 
system  was  far  removed  from  equilibrium  solute  distribution  conditions. 
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